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Abstract. Wo adapt a simplified version of the Multi-Scale Analysis presented 
in [8] to multi-particle tight-binding Anderson models. Combined with a recent 
eigenvalue concentration bound for multi-particle systems [7], the new method 
leads to a simpler proof of the multi-particle dynamical localization with more 
optimal decay bounds on eigcnfunctions than in I4ll5lll0l , for a large class of strongly 
mixing random potentials. All earlier results required the random potential to be 
IID. We also extend the result on multi-particle localization to models with a 
rapidly decaying interaction. 

Note: this is an improved version of the manuscript originally up- 
loaded on 11.06.2011. 



1. Introduction. The model, assumptions and the main results 

Analysis of localization phenomena in multi-particle quantum systems with nontriv- 
ial interaction in a random environment is a relatively new direction in the Anderson 
localization theory, where during half a century, since the seminal paper by P. Anderson 
PQ , most efforts were concentrated on the study of disordered systems in single-particle 
approximation, i.e., without inter-particle interaction. While in numerous physical 
models such an approximation is fairly reasonable, it was pointed out already in the 
first works by Anderson that multi-particle models presented a real challenge. 

The number of results on multi-particle localization obtained both by physicists and 
mathematicians remains yet rather limited. We do not review here results obtained by 
physicists, based on methods of theoretical physics and considered as firmly established 
by the physical community. The rigorous mathematical results on multi-particle local- 
ization obtained so far (cf. [5hTT|. [U[S]) apply to TV-particle systems with arbitrary, 
but fixed N > 1, and the range of parameters (such as the amplitude of the disorder 
and/or proximity to the edge(s) of the spectrum) is rapidly degrading as N — > oo. 

In this paper, we study spectral properties of random lattice Schrodingcr operators 
(LSO) in the general framework of the Multi-Scale Analysis (MSA) developed in [T3l 
rHlfTBlrM] . Specifically, we study random Hamiltonians of the form 

H(w) = H + gV(w) + U, g g K, (1.1) 

where Ho is a finite-difference operator representing the kinetic energy, e.g., the nearest- 
neighbor lattice (negative) Laplacian (—A), U is the operator of multiplication by the 
interaction potential x M> U(x), and V(w) is the operator of multiplication by the 
function x M> V(x; w) = V{x\;oS) + ■ ■ ■ V(x^;cj), where V : Z d x O — > K is a random 
field relative to some probability space (17, J 7 , P). As was shown in [TU], [U[S], Anderson 
localization can be established in the entire space 7-Ln, but in this paper we consider 
only the restriction of H(w) to the fcrmionic subspacc. The bosonic subspace can be 
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treated essentially in the same way. Using either kind of quantum statistics makes the 
scaling analysis more transparent and proofs substantially simpler. 

Aizenman and Warzel [H[S] established the exponential dynamical localization in 
multi-particle systems in the framework of a parametric analysis allowing to prove local 
stability (in the parameter space) of the Anderson localization phenomenon (including 
exponential strong dynamical localization) for Anderson-type Hamiltonians with an 
IID external random potential under perturbations by a short-range interaction. 

Apart from simplifications of the Multi-Particle Multi-Scale Analysis (MPMSA) 
developed in [9j[T0], the novelty of this paper is three- fold: 

• the external random potential is not assumed to be IID, but is strongly mixing; 

• the interaction potential is not necessarily of finite range, but can be exponen- 
tially or even sub-exponentially fast decaying at infinity; 

• the decay of eigenfunctions is proven with respect to a norm-distance in Z Nd 
(more precisely, in a symmetrized norm-distance), while earlier results used 
explicitly ([4]) or implicitly (flOp a more complicated and less natural notion 
of distance in the multi-particle configuration spac^H- 

Note that our proof of multi-particle localization for strongly mixing random po- 
tentials and rapidly decaying infinite-range interaction can be extended to the entire 
Hilbert space "Hat, at the price of greater complexity On the other hand, the proofs can 
be made simpler and more streamlined, if only a power-law decay of the eigenfunction 
correlators were to be established (as is customary in the framework of the MSA); the 
latter is not to be confused with the decay of eigenfunctions which is exponential. 

We consider here only strongly disordered systems; an adaptation of our approach 
to localization at "extreme" energies in weakly disordered multi-particle systems, as 
well as to weak perturbations of localized non-interacting multi-particle systems, has 
been recently obtained by Ekanga |17j in the framework of his PhD project. 

The new scaling procedure used in this paper is an adaptation of the method pro- 
posed earlier in [5] in the single-particle context; it is close in spirit to a very simple 
method developed by Spencer [24] for a fixed-energy analysis of Green function:^. 

The proofs of all statements not given in the main text can be found in Appendix. 

The structure of this paper is as follows. 

• In Section [5] we describe main geometrical notions and constructions relative 
to fermionic systems; for brevity of presentation, we consider first the case 
of one-dimensional particles. Nevertheless, our scheme in Sections [3H5] is not 
specific to one dimension; indeed, it has to be stressed that from the analytic 
point of view, an A^-particle system in Z 1 (starting already with N = 2) gives 
rise to a multi-dimensional (discrete) Schrodinger operator. 

• In Section [3] we describe the analytic core of the new scaling procedure, which 
is essentially the same as in the single-particle case treated in [3]. 

• The probabilistic inductive bounds for the model with a finite-range interaction 
are given in Section [4] The key statements here are Lemma 14.41 Lemma 14.51 
and Theorem 14.21 A reader familiar with |10) may notice that an important 
component of the the proof of Lcmma [4.4[ given by Lemma l4.31 is much simpler 



A detailed analysis of the nature of this technical difficulty, performed by Aizenman and Warzel 
[4], proved instrumental for a solution proposed later in [7], 

2 I thank Tom Spencer, Abel Klein and Boris Shapiro for fruitful discussions of the work 1241 . 
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than its counterpart in |10j . due to the use of eigenfunctions instead of Green 
functions. 

• An adaptation to infinite-range interactions is given in Section [5j 

• The derivation of the strong dynamical localization from the results of the 
scaling analysis is given in Section [51 where we follow the same path as in [5] . 
A reader familiar, e.g., with [T51[5D] may want to skip Section [5] 

• An adaptation to more general graphs is discussed in Section [7J 

• For the reader's convenience, we prove in Appendix the new eigenvalue con- 
centration (EVC) bound (Theorem 14.11 cf. also [7]), allowing to establish the 
exponential decay of eigenfunctions with respect to the max-norm in the N- 
particlc configuration space. 

To conclude the introduction, note that we do not discuss here the ergodicity issues 
for the Hamiltonians H(w), for several reasons. 

First of all, it is well-known already in the single-particle localization theory, where 
the Hamiltonian has the form H{uj) = Ho+gV(uj), say, on the lattice Z d , the ergodicity 
of the random field V : 7L d x f2 — > R per se is irrelevant for the usual proofs of complete 
localization (pure point spectrum) under the assumption of strong disorder (\g\ ^> 1). 
In the case of an IID potential, the location of the a.s. spectrum can be described 
withthe help of the Weyl argument, again, without using explicitly the ergodicity (which 
is, of course, granted for an IID random field). 

The second reason is that already in the case of a two-particle Hamiltonian H(w) = 
H + W(x; oj) = H + V(x; w) + U(x), the random field 

((xi,x 2 ),oj) H> V(xi;oj) + V(xi;lj) + U(xi,x 2 ) 

is no longer stationary. Specifically, let d = 1 and V(x;uj) € [0,1], say, uniformly 
distributed in [0,1], and U(a;i,a;2) = U l{ Xl =x 2 }> U < —2. Then for all points of 
the form x = (x,x), W(x;w) < 0, while for all others W(x;w) > 0. The translation 
invariance of the field W(x; lj) holds only for a subrgoup of "diagonal" shifts {x\, x 2 ) >— >• 
{x\ + a, X2 + a), a € Z. 

The third, more important and less formal reason is that the language of the density 
of states (DoS), very helpful and instructive in the context of single-particle models, is 
much less so in the framework of multi-particle models with a nontrivial interaction. 
Recall that, according to an earlier result by Klopp and Zenk [33], proven for multi- 
particle Schrodingcr operators in R d with decaying interaction, the DoS is the same 
for operators with and without interaction. (An adaptation of their techniques to 
the lattice is not difficult.) Consider again the two-particle operator H(w) from the 
previous paragraph; let Ho > 0. Then the density of states is supported by the non- 
negative half-line M. + , but it is obvious that , with U < — | Ho 1 1 — 1 1 V 1 1 , the spectrum has a 
negative component, for the quadratic form (/|H/) is not sign-definite, having positive 
and negative diagonal matrix elements in the delta-basis. Speaking informally, the DoS 
measure indicates the location of the "bulk" spectrum, while the negative interaction 
gives rise to some "internal surface" spectrum due to eigenfunctions (square-summablc 
or generalized) essentially supported by a neighborhood of supp U and decaying away 
from it (which can be seen, e.g., with the help of the Combes-Thomas argument). 

The sign of the interaction is not crucial for the above observation; this can be seen 
already from the fact that finite-difference Schrodingcr operators ±H(w) have similar 
qualitative spectral properties. Specifically, let U(a;i,a;2) = Ul[ Xl=X2 y, U = 11 = 
1 + ||H || + ||V||, H > 0, ||H || = 8. Then the Dos is the same as with U = 0, hence, 
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contained in [0, ||H || + ||V|| ] = [0, 10] (recall that < V < 1 in this example). On 
the other hand, for x = (0, 1), one has, by non-negativity of Ho + V, 

<1 X |H1 X ) > 17 = 11. 

Therefore, the spectrum of H(w) in (10,11] is nonempty, while the DoS vanishes in 
(10,11]. 

Nevertheless, it is true that, even though the potential random field W(x;oj) is 
not ergodic, Hamiltonians H(oj) still form an ergodic family of operators, hence their 
spectral components are a.s. nonrandom. 

2. Basic definitions, assumptions and main results 

2.1. Configurations of indistinguishable particle in Z 1 . In the first part of this 
paper, wc work with configurations of N > 1 quantum particles in the one-dimensional 
lattice Z. In quantum mechanics, particles of the same kind are considered indistin- 
guishable; more precisely, depending on the nature of the particles, the wave functions 
describing N > 1 particles have to be either symmetric (Bose-Einstcin quantum sta- 
tistics) or antisymmetric (Fermi-Dirac quantum statistics). We choose here fermionc 
case; this gives rise to slightly simpler notations and constructions. 

For clarity, we use boldface symbols for objects related to multi-particle systems. 

Quantum states of an iV-particle fermionic system in Z are elements of the Hilbcrt 
space H N = H N '~ of all square-summable antisymmetric functions ^ : Z, N — > C, 
with the inner product (■ | •) inherited from the Hilbert space £ 2 (Z N ) = (H 1 )® 1 ^ . In 
this particular case where the "physical" configuration space is one-dimensional, H 
admits a simple representation which we will use. 

First, note that any antisymmetric function x = (x\, ■ ■ ■ , Xjv) >-> *&{xi, . . . , Xn) 
vanishes on all hypcrplanes {x : a;, = Xj}, 1 < i < j < N. Further, any function 
defined on the "positive sector" Z™ = {(x\, . . . ,xn) ■ x\ > X2 > • ■ ■ > xn} admits a 
unique antisymmetric continuation to Z := 1 N . An ortho normal basis in ft can 
be chosen in the usual form 

* a (x) = -L V (-iy ® l a ae2" #{a u ...,a N } = N. 

Here &n is the symmetric group acting in 7L N by permutations of the particle posi- 
tions, 7r(a) = (a^-i(i), . . . , a 7r - 1 (iv))) an( l ( — l) w £ 1} denotes the parity of the 
permutation tt. It is readily seen that H N is unitarily isomorphic to the Hilbert space 
( 2 _{Z™) of square-summable functions on Z™; equivalcntly, one can consider square- 
summable functions on the set Z™ = {(xi > ■ • • > xjy} vanishing on the boundary 
Z™ = {x € Z™ : 3 i ^ j, Xj = Xj}. Indeed, the isomorphism is induced by the 
bijection between the orthonormal bases {<& a } and {l a }: 

* a <-> l a = l ai ® • l a „, a€Z™. 

The subspace H N is invariant under any operator commuting with the action of the 
symmetric group ©at. 

An advantage of the above representation is that Z™ inherits its explicit, natural 
graph structure from Z N . 



DIRECT SCALING ANALYSIS. II. MULTI-PARTICLE SYSTEMS 



5 



Occasionally we will denote by x the ordered configuration (i.e., a vector) {x\, . . . ,xn) 
corresponding to an unordered configuration x = {x\, . . . , xn}. 

2.2. Fermionic Laplacians. Any unordered finite or countable connected graph (Z, 8) 
with the set of vertices Z and the set of edges 8 is endowed with the canonical graph 
distance (x,y) i-> dz(x,y) (defined as the length f the shortest path x ~+ y over the 
edges) and with the canonical (negative) graph Laplacian (— Az)' 

{-A z f){x) = £ (/(*) - f(y)) = n z (x)f(x) £ f(y); 

(x,y) {x,y) 

here (x,y) denotes a pair of nearest neighbors, and nz(x) is the coordination number 
of the point x m Z, i.e., the number of its nearest neighbors. 

In particular, one can take Z = 7L N with the edges (x, y) formed by vertices x, y 
with | a; — y\i := \x± — yi\ + ■ ■ ■ + \xn — J/jv| = 1- I n other words, the vector norm | • |i 
induces the graph distance on Z . The (negative) Laplacian on Z w , which we will now 
denote by (—A), can be written as follows: 




where J.W is the identity operator acting on the i-th variable, and (— A^) is the 
one-dimensional negative lattice Laplacian in the j-th variable: 

(-A«/)fe) - 2f( Xj ) - f( Xj - 1) - f( Xj + 1), Xj g Z. 

Its restriction to the subspace H N of antisymmetric functions can be equivalently 
defined in terms of functions supported by the positive sector Z™ , hence vanishing on 
its border Z™ . Indeed, the matrix elements of A in the basis 3? a can be nonzero only 
for pairs <& a , 3>b with |a — b|i = 1, so that, for some j £ [1, N], \ctj —bj\ = 1, while for 
all i ^ j, ai = In. If a, b e Z», then (* a | A | * b ) = ^<l a | A | l b ). If, say, a € Z", 
then the respective matrix element of the Laplacian's restriction to Z™ with Dirichlct 
boundary conditions on Z™ vanishes, but so does the function 4> a (which is no longer 
an element of the basis in H N ). Therefore, up to a constant factor, the restriction of 
the TV-particle Laplacian to the fermionic subspace H is unitarily equivalent to the 
standard graph Laplacian on Z™ . From this point on, we will work with the latter, 
occasionally making references to the space t 2 {X N ). 

It will be convenient to use in the course of the scaling analysis a different notion of 
distance on 7* N (hence, also on Z™ C Z N ): the max-distance defined by 

P(^y) = max ■ d z (xj,yj), 

1<J<N 

and work with balls relative to the max-distance, 

N 

Bl(x) - {y : p(x,y) < L} = X B L ( Xj ). (2.1) 

3=1 

Here B L (x) = [x - L, x + L] n Z. 

In the positive sector Z™ , the above factorization of p-balls is subject to the condi- 
tion that the RHS of (|2.1j) is itself a subset of the positive sector Z™ (but the inclusion 
"C" always holds true). Considering configurations x = {x\, . . . ,x n } as subsets of Z, 
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one can define the distance between two configurations x' G -H> and x" g ^ n,,) in 
a usual way: 

dist(x',x") = min min \u— v\. 

Lemma 2.1. Lei x G Z™ be a union of two subconftguraitons x' G Z^ 1 ^ and x" G 
.Z^ 1 \ such that dist(x',x") > 2L. Then the following identity holds true: 

Bf ) (x)=Bi"')(x')xBi"" ) (x"). 
The proof is straightforward and will be omitted. 

The graph distance d Z N will be useful in some geometrical constructions and defi- 
nitions, referring to the graph structure inherited from Z N . 

Given a subgraph A C Z™ , we define its internal, external and the so-called graph 
(or edge) boundary, in terms of the canonical graph distance: 

d-A={yeA: d^(y, A c ) = 1 }, d+A = d~ A c 
OA={(x,y)eAxA c : d^(x,y) = l}. 

We also define the occupation numbers relative to a configuration x = (x\, . . . ,xjy). 
Namely, define a function n x : Z H> N by n x (y) = #{j : Xj = y}, y G X . 

2.3. Multi-particle Hamiltonians. The matrix elements of resolvents G\{E) = 
(Ha — E)~ x , for E S(Ha) = spcc(Ha), in the canonical delta-basis, usually referred 
to as the Green functions, will be denoted by Ga(x, y;E). In the context of ran- 
dom operators, the dependence upon the element u> G fl will be often omitted, unless 
required or instructive. Similar notations will be used for infinite A. 

2.4. Graph Laplacians and fermionic Hamiltonians. Given a real-valued func- 
tion W : Z™ M, we identify it with the operator of multiplication by W and 
consider the fermionic iV-particle random Hamiltonian 

H(w) = H + W(w), 

where Ho is a second-order finite-difference operator in £ 2 (Z™), for example the graph 
Laplacian A. We assume that 

W(x) = W(x; u) = ffV(x; u) + U(x) 

where g~V(-;u)) is the external random potential energy of the form 

V(x; u) = V(xi; u) + ■ ■ ■ + V(x n ;uj), 

and V : Z x — >• M is a random field on Z relative to a probability space (f2,#, P); the 
expectation relative to the measure P will be denoted by E [ ■ ] . Our assumptions on V 
arc listed below (cf. Wl- W3). 

Further, U is the interaction energy operator; for notational brevity, we assume that 
it is generated by a two-body interaction potential U, 

U(x)=J2UQx i -x j \), 

satisfying one of the hypotheses U0-U1 (see below). 

Note that, without loss of generality, if ||Ho|| < oo, then one can always assume that 
the absolute values of the matrix elements of Hq are bounded by 1. Otherwise one can 
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consider a two-parameter family Ho + /i _1 (jV(x; uS) + U(x)) = h^ 1 (hUo + <?V(x; uS) + 
U(x)), so the two operators share the eigenvectors, while the matrix elements of ft-Ho 
can be made smaller than 1 by taking h small enough. 

One can also consider higher-order finite-difference operators Ho; this requires only 
minor technical modifications; see the discussion in Section [3] 

2.5. Assumptions on the random potential. We assume that the random field 
V : Z d x fi — > R is (possibly) correlated, but strongly mixing; this includes of course 
the IID potentials treated earlier in |51ITU] and gJ[S]. Lct F v . x (t) = ¥{V(x;lu) <t}, 
Z d , be the marginal probability distribution functions (PDF) of the random field V, 
and Fv, x (t \ d^x) '■= P { V(x; uj) < t \ } the conditional distribution functions (CDF) 
of the random field V given the sigma-algebra 3s^ x generated by random variables 
{V (y; ui),y 7^ x}. Our assumptions on correlated potentials are summarized as follows: 
Wl The marginal CDFs are uniformly Holder-continuous: for some k > 0, 

esssup sup sup[F v , x (a + s\^jt x ) - F v , x (a\$^ x )\ < Const s K . (2.2) 

W2 (Rosenblatt mixing) For any pair of subsets B', B" C 1> d with d(B',B") > L, any 
events £' € #b' , £" <= 5b" and some C > 

{£' C\£"} -¥{£'}¥{£"} <e- CXv?L . (2.3) 

One can easily check that for any p > 0, a € (1, 2), b G (0, 1), arbitrarily large a > 
and sufficiently large Lq £ N 

e -Cln 2 L° < / £ a*} ^ ^.^Q (2 4) 



The rate of decay of correlations indicated in the RHS of (|2.4p is required to prove 
dynamical localization bounds with decay rate of EF correlators faster than polynomial, 
and it can be relaxed to a power-law decay, if one aims to prove only a power-law decay 
of EF correlator^. 

Assumptions W1-W2 arc sufficient for the proof of spectral and strong dynamical 
localization in multi-particle systems. In particular, the role of Wl is to guarantee 
Wcgner-type estimates used in the MPMSA scheme. However, it was discovered in [?J[5] 
and in [SKTU] that traditional, Wegner-type EVC estimates do not provide all necessary 
information for efficient decay bounds on the eigenfunctions of multi-particle operators. 
More precisely, conventional EVC bounds seem so far insufficient for the proof of the 
exponential decay of eigenfunctions with respect to a norm in the configuration space 
of TV-particle systems, starting frorrQ N = 3. For this reason, we proposed earlier [7] 
a new method allowing to compare spectra of two strongly correlated multi-particle 
subsystems and extending the Wegner-type EVC estimate. For the new method to 
apply, one needs an additional assumption on the random potential field which we will 
describe now. 

First, introduce the following notations. Given a lattice parallelepiped Q C Z d , we 
denote by £q(w) the sample mean of the random field V over the Q, 

tQ(u) = (V) Q = \Q\- 1 J2VM 
xeQ 



In this case, Wl can be relaxed to a form of log-Holder continuity of the marginal CDFs. 
^For N = 2, one can use a simpler EVC bound proven in [9]. 



s 



V. CHULAEVSKY 



and define the "fluctuations" of V relative to the sample mean, r\ x = V(x,cj) — 
£q(lj), x £ Q. Denote by dv.Q the sigma-algebra generated by {rj x , V y : x G Q,y G" Q}, 
and by F^(- |3V,q) the conditional distribution function of £q given 3v,Q- Further, 
consider parallelepipeds Q C Z d with diam(Q) < R < oo, and introduce an $v,Q- 
measurable random variable 

v r (s;uj) := sup\F s (t + s\$ VtQ ) - Ft(t\$ VtQ )\, s > 0. (2.5) 
teR 

We will assume that the random field V fulfills the following condition: 
W3: There exist C , C" , A', A", b' , b" G (0, +oo) such that 

Vse[0,l] ess sup P { v R (s; u) > C'R A 's b ' }< C"R A "s b ". (2.6) 

In the particular case of a Gaussian IID field, e.g., with zero mean and unit variance, 
£q is a Gaussian random variable with variance \Q\ 1 independent of the "fluctuations" 
rj x , so that its probability density JVjqi is bounded, although ||2V,|Q| ||oo ~ IQI 1 ^ 2 — > oo 
as |Q| — > oo. The property W3 has been recently proven for a larger class of IID 
potentials by Gaume [T^] in the framework of his PhD project. Note that the property 
W3 is not quite obvious, since the conditional distribution given by F^(t\^Sv,Q) can 
be very singular - even discontinuous - for some conditions, e.g., in the case of an IID 
random field V with uniform marginal distribution. Nevertheless, the proof of W3 for 
an IID random field with uniform marginal probability distribution is quite simple. 

2.6. Assumptions on the interaction potential. We assume that the interaction 
potential U generating the interaction energy U satisfies one of the following decay 
conditions: 

U0 There exists ro < oo such that 

Vr>r U(r)=0. (2.7) 
Ul There are some C G (0, +oo), S G (0, j^^j and 9 G (0, j^g) such that 

Vr>0 \U(r)\ < Ce~ crl ~ e . (2.8) 

Naturally, U0 =>■ Ul. Wc will prove the multi-particle localization first under the 
strongest assumption U0 (leading to a simpler proof), to illustrate the general structure 
of the new MPMSA scheme, and then extend the proof to interactions satisfying Ul. 

2.7. Main results. 

Theorem 2.1. Assume that the random field V fulfills conditions W1-W3, and the 
interaction U fulfills one of the conditions U0, Ul. There exists go G (0, +oo) such that 
if \g\ ^ <?0; then with probability one, the fermionic operator H(w) = A + gV(cj) + U 
has pure point spectrum and all its eigenf unctions ^fj(uj) are exponentially decaying at 
infinity: for each tyj, some and all x with p(x,x) large enough, 

|* 3 -(x,w)| < e -"v(*.*), m > 0. (2.9) 

In addition, for all points x, y with p(x, y) large enough and some a,c > 0, for any 
bounded Borel function f : M. — > M 

E[|(l y |/(HH)| l x )|] ^e-^^WfU. (2.10) 
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Consequently, for any finite subset K C Z 



Nil 



E 



e aln 



l + c 



X /(HH)1 K 1 <C(K)||/|| 00 <oo (2.11) 



where operator Xk is £/ie defined by (Xr/)(x) = (p(K,x) + l)/(x). /n particular, 



E 



sup 



e e v ' Ik 



< oo. (2.12) 



3. Deterministic bounds 

3.1. Geometrical resolvent inequality. The most essential part of the scaling anal- 
ysis concerns finite-volume approximations H\(u)) of the random Hamiltonian H(w), 
acting in finite-dimensional spaces H\ := £ 2 (A), A C Z , |A| = card A < oo. 
Operator A can be represented as follows: 

-A = n zN - ( r xy + r yx ), (r xy /)(x) := 5xy/(y), 
(x,y> 

and S xy is the Kronecker symbol. Similar formulae are valid for the restriction Aa of 
A to a finite subset A C Z N with Dirichlct boundary conditions outside A; in this 
case, one has to keep only the pairs (x, y) € A x A. 

We denote by G\(E) = (H A - E^ 1 the resolvent of H A and by G(x,y;E) the 
matrix elements thereof (a.k.a. Green functions) in the standard delta-basis. 

The so-called Geometric Resolvent Inequality for the Green functions can be easily 
deduced from the second resolvent identity: 

|G BL (x,y;£)| <Q max |G Bf(x) (x, v;£)| max |G B[ (V, y; E)\ (3.1) 
vea-B ( (x) v'ed+B f (x) 

where 

C t = \dB t (v)\ < C(N, d)£ Nd . (3.2) 
Clearly, (|3.1[) implies the inequality 



\G BL ( X ,y;E)\<{C ( max |G Bf(x) (x,v;£)| ( max |G Bi (v'j;£)| (3.3) 

which is weaker than (|3.1[) . but sufficient for the purposes of the scaling analysis. 
Sometimes we shall use the following inequality which stems from (|3.ip : 



|G Bi (x,y;£)| <C,||G B , (u) (i;)|| max |G Bi (v, y; E)\ . (3.4) 

v:p(u,v)<«+l 

Similarly, for the solutions ip of the eigenfunction equation Hi/? = Eif) we have 

|^(x)| < Q ||G B , (x) (£)|| max |t/>(y)|, (3.5) 

y:p(x,y)<«+l 

provided that E is not an eigenvalue of the operator H Bf ( u ) and x € B^(u). 

3.2. Radial descent bounds for Green functions and eigenfunctions. The an- 
alytic statements of this subsection apply indifferently to any LSO, regardless of their 
single- or multi-particle structure of the potential energy operator. To avoid any confu- 
sion, we denote the dimension of the lattice by D; in applications to A^-particlc models 
in Z d , one has to set D = Nd. 
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Definition 3.1. fix an integer £ > and a number q G (0,1). Consider a finite 
connected subgraph A c IP and a ball Br{u) C A. A function / : Br(u) — > R+ is 
called (£, g)-subharmonic in Br(u) if for any ball Bg(x) C Ar(u) one has 

f(x) < q max f{y). (3.6) 

y:d(a;,y)<« 

Lemma 3.1. Le be given integers L > £ > and a number q € (0,1). If f : A — > M+ 
on a finite connected graph A is (£,q)-subharmonic in a ball Bl(x) C A, i/ien 

/(x)<gL^J^(/,B)<^M(/,B), M(/,A):=max|/(a:)|. (3.7) 

Proof. The claim follows from (|3.6[) by induction ( "radial descent" ) : 

f(x) < qM(f,B t+1 ) <■■■< q j M(f,B j( z +1) ) <■■■< ql^\M(f,B). □ 

Lemma 3.2. Let A be a finite connected graph and / : A x A -) IR + . (x'.x") t— > 
f(x',x"), a function which is separately (£,q)-subharmonic in x' G B r i(u') C A and in 
x" G B r » (u") C A, with d(u', u") > r' + r" + 2. T/ien 

f(u',u") < gL^J + l.T#l<M(/,B) < q r ' + '-+i 2e M(f 7 B). (3.8) 
Proof. Fix any point y" G B r " + i(u") and consider the function 

f y '.:y'^f{y',y") 
which is {£, g)-subharmonic in B r i(u') C A, so by Lemma |3. 11 

/»»(«') = /(«',!/") < J^M/,Ax A). (3.9) 

Now introduce the function /„< : y" M> f(u', y") which is (£, <7)-subharmonic in B r » (u") C 
A and bounded by the RHS of (|3.9|) . Applying again Lemma |3~T1 the claim follows. □ 
The relevance of the above notions and results is illustrated by the following 

Lemma 3.3 (Cf. Lemma 3.3 in [5]). Consider a ball Bl(u) and an operator H = 
Hb l (u) with fixed (non-random) potential V. Let {ipj, j = 1, . . . , |B/,(u)|} be the nor- 
malized eigen) 'unctions of H . Pick a pair of points x' , y' G Bl(u) with d(x' , y') > 2(£+l) 
and an integer R G [£+2, d(x' ,y') — (£+2)]. Suppose that any ballB((v) withv G B^(x') 
is (£7, m)-NS, and set 

q = q{D,£;E) = C(D)e-^ m ^ e (3.10) 
with the constant C(D) defined in the same way as C{d) in (|3.2[) . Then: 

(A) the kernel Hip.(x,y) of the spectral projection Li^ . = \tpj)(ipj\ is (£+1, q)-subharmonic 
in x G Br(x'), with global maximum < 1; 

(B) i/Bi(«) is also _E-NR, then the Green functions G(x,y;E) = (H — £') _1 (x, y) are 
{£ + 1, q)-subharmonic in x G Bfl(x'), wii/i global maximum < e L . 

Remark 3.1. The general idea of "two-sided" bounds on the Green functions (as func- 
tions of two arguments) can be found already in the Spencer's paper [24] . This allows 
to avoid a more complex geometrical construction going back to [T3], required when 
more than one "unwanted", "singular" locus is to be allowed in the course of the scaling 
step. 



Remark 3.2. An extension of Lemma 13.31 to finite-difference kinetic energy operators 
Ho of higher order is quite straightforward. Specifically, if Ho has range £o > 1, i.e., 
(l x | H l y ) = for p(x, y) > £ , the denominators £ + 1 in Eqns (|3T7|) . (pO]) . ([33)) arc 
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to be replaced by i + £q. Recall that we assume (implicitly) that the matric elements 
of Ho are bounded by 1 (as they are for the graph Laplacian), but the reduction to 
this case can be made by rescaling the operator Ho, i.e., by replacing Ho hH-o, with 
h > small enough (provided that ||Ho|| < oo). Apart from a modification of some 
auxiliary constants, this has no impact on the qualitative final result on multi-particle 
localization under the assumption of strong disorder. 

3.3. Localization and tunneling in finite balls. 

Definition 3.2. Given a sample V{-;uS), a ball B^(u) is called 

• E'-non-resonant (£?-NR ) if ||Gb(-E; cj)j| < e +L ' 3 , and i?-resonant otherwise; 

• completely E- non-resonant (S-CNR) if it does not contain any E-K ball B^(u) C 
Bi(u) with £ > L 1 / 01 , and -E-partially resonant (E-PK), otherwise. 

Definition 3.3. Given a sample V(-;u>), a ball Bj,(u) is called (E, m)-non-singular 
((£?,m)-NS) if 

max \G(u iy ;E-iu)\<e-^ m ' L)L+2L ^ (3.11) 
v edB L (u) 

where 

7(m,L) :=m(l + L- T ) , r=l/8. (3.12) 
Otherwise, it is called (E, m)-singular ((E,m)-S). 

In the next definition, we use a parameter g = (a — l)/2; with a = 4/3, one obtains 
£i=l/6 and (1 + g)/a = 7/8. 

Definition 3.4. Given a sample V(-;oj), a ball Bj,(u) is called m-localized (m—Coc, 
in short) if any eigenfunction ipj of the operator Hb l ( u ) operator satisfies 

|^.(x)^(y)| <e^ WWw » (3.13) 

for any pair of points x, y € B^(u) with p(x, y) > = L 7 ^ s . 

Definition 3.5. We will call a pair of balls B^(x), B^(y) distant iff p(x, y) > CatL, 
with Cat = 2A N + 3, A N = 4N. 



(The role of the constants Am and Cn will become clear in subsection I4.3H 

Definition 3.6. A ball B£ Q (u) is called m-tunneling (m-T) if it contains a pair of 
distant m—MCoc balls B^(x), B^(y), and m-non-tunneling (m-NT), otherwise. 

It is to be emphasized that, unlike the property of .E-resonance or (E, m)-singularity, 
the tunneling property is not related to a specific value of energy E, and tunneling 
even in a single ball occurs with a small probability. 

Lemma 3.4. Let E E M. and an m—Coc ball Bj,(u) be given. //Bi(u) is also £7-NR, 
then it is (£',m)-NS. 



Proof. The matrix elements of the resolvent Gb (E) can be assessed as follows 

(x)IIV> 

\E-E. 



G B (x,y;E)|< 2^ \E-E\ ' ^ ^ 



If dist(£7, ct(H-b)) > e l>3 and ln(2Z + l) d < L 13 , then the m—Coc property implies 

|G B (x,y;£)| < e ^(^L)L+L^+ln\-B\ < e - 7 (m,L)L+2L^ Q 
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Observe that with m > 1 and 1 — r > /?, for Lq large enough 

m(l + r )L fc - 2if = mL k + mL k ~ T - 2Ll > m(l + \L k T )L k . (3.15) 

We will see that the condition m > 1 is fulfilled for |g| large enough; cf. Lemma \4. II 

From this point on, we will work with a sequence of "scales" - positive integers 
{Lk, k > 0} defined recursively by L k+ i = \L\V\ , L > 2. For clarity, we keep the value 
a = 4/3; observe that a 2 < 2. In several arguments we will require the initial scale Lo 
to be large enough. 

Lemma 3.5. There is L (1) G N such that for L > 

(A) if a ball B£ fc+1 (u) is .E-CNR and contains no pair of distant (E, m)-S balls 
Bz, fc (v), Bi fc (w), then it is also (E,m)-NS; 

(B) if a ball B Lfc+1 (u) ism-NT and E-CNR, then it is also (£,m)-NS. 

Proof. (A) By assumption, cither B/, fc+1 (u) is E'-CNR and contains no (E, m)-S ball 
of radius L k , or there is a point w € Bi fc+1 (u) such that any ball Bi t (v) C Bi t (u) 
with p(w, v) > C N L k is (E, m)-N^. 

In the former case, such an exclusion is unnecessary, but in order to treat both 
situations with one argument, we can formally set w = u (or any other point). 

Fix points x, y with R :— p(x, y) > L k +Q . By triangle inequality, 

p(x, B {CN „ 1]Lk (w)) + p(y, B {CN _ 1)Lk (w)) > R - (2C N - 2)L k . 
Assume first that 

r' :=p(x,B (Cw _ 1)ifc (w)) >L k + l, 

r" :=p(y,B (CN _ 1)Lk (w))>L k + l. 

All balls of radius L k both in B r / (x) and in B r » (y) are automatically disjoint from 
B2L fc (w), thus (E, m)-NS. Furthermore, 

r> +r" >R- 2(C N - l)L k -2>R- 2C N L k . 

Consider the set B = B r /(x) x B r "(y) and the function / : B — > C defined by 
/(x', x") = Gb L;s+1 (x', x"; E). Since E is not a pole of the resolvent Gb £ , Jj+i (u)(')i ^ IS 
well-defined (hence, bounded, on a finite set). By Lemma l3~3l / is (L k , q)-subharmonic 
both in x' and in x", with q < e -7(m,L h ,n) > Therefore, one can write, with the conven- 
tion — InO = +oo, using Lemma 13.21 and setting for brevity K := M — n + 1: 

R — 2C N L k — 2L k . 



-ln/(u,y) > - In 



e k + 1 



= mR 



l & +i^V x o^jv^ ^ k ) ^ k+1 ^3 16 ^ 



(i + i^r) (i - (i - ZC N L h 



mR 



> m (1 + |Lfc r ) # > lK ife+i)^- 
If r' = (resp., r" = 0), the required bound follows from the subharmonicity of the 
function J(x',x") in x" (resp., in x'). 

(B) Assume otherwise. Then, by assertion (A), the _E-CNR ball Bt t+1 (u) must 
contain a pair of distant (E, m)-S balls Bi t (x), B^^y). Both of them are E-NR, 



^The value of CV is stipulated in Definition 13.51 but it is irrelevant for this proof, provided that 
Lq is large enough. 
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since Bi t+1 (u) is .E-CNR. By virtue of Lemma 13.41 both Bi t (x) and Bi fc (y) must be 
m—NCoc, so that Bt t+1 (u) must be m-T, which contradicts the hypothesis. □ 

Lemma 3.6. There is G N such that, for all L > IA 2 ), if for any E G R a ball 
Bi t+1 (u) contains no pair of distant (E,m)-S balls of radius Lk, then it is m—Coc. 

Proof. One can proceed as in the proof of the previous lemma, but with the functions 
fj : (v',v") H> \ipj(v')ipj(v")\, where t/> ■ , j G [1, |Bi fc+1 (u)|], are normalized eigenfunc- 
tions of operator HB lfc i ( u ). Notice that the i?j -non-resonance condition for B^ fc+1 (u) 
is not required here, since \\ipj\\ = 1, so the function \fj \ is globally bounded by 1. Let 
x, y G Bi k+1 (u) and assume that p(x, y) = R > L^ +e . 

Arguing as in the previous lemma and using the subharmonicity of the function 
fj{- , ■) (cf. Lemma l3~3"|) . we obtain, for any pair of points with p(y',y") > L\ +e : 



ln|^.(x)^(y)| > - In 



^ e -m(l + ^L-^)L k y 



(3.17) 



A direct comparison with the RHS of the first equation in (|5.4[) shows that the RHS of 
(I3.17[) is bigger, owing to the absence of the factor e L ^ . Therefore, it admits the same 
(or better) lower bound as in in (|5.4[) . □ 
Below wc summarize the relations between the parameters used in the above proof: 

p(x,y) >Ll±l, 0<T<p,l + g<a;(3<l-T. (3.18) 
These conditions are fulfilled, e.g., for 

a = 4/3, ,0 = 1/6, r = l/8, [3 = 1/2. (3.19) 

4. Scaling analysis of eigenfunctions. Finite- range interaction 

4.1. Initial scale estimates. 

Lemma 4.1 (Initial scale bound). For any Lq > 2, m > 1 and p > there is go = 
go(m,p, Lq) < oo such that for all g with \g\ > go, any ball Bi (u) and any £eM 

P{B io (u) is(E,m)-S} <L~ P (4.1) 

P{B io (u) is m-MCoc] <L~ P . (4.2) 

See the proof in Appendix. 

4.2. Eigenvalue concentration bounds for distant pairs of balls. 

Theorem 4.1 (Cf. [7])- Let V : Z x O — !• R be a random field satisfying W3. For all 
L > large enough, any pair of balls B^(u'), B^(u") with p$(u',u") > ANL the 
following bound holds for any s G (0, 1]: 

P{dist(a(H Bl ( u0 ),a(HB l( u")))<s} < (2L + l) 2Nd h L (2s) (4.3) 

where 

h L (s) := C"L A "(2sf + |B L (u')| |B L (u")| C'L A '(2sf. (4.4) 
Consequently, this implies the following bound for ft 1 G (0, f3), all L large enough and 
any pair u' , u" , with d(u',u") > ANL: 

P{3£gR: B L (u') and B L (u") are E-PR} <e- L ^' (4.5) 
Sec the proof in Appendix. 
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4.3. "Partially interactive" A-particle balls. Unlike the situation in single-particle 
models where the configuration space (e.g., the lattice Z d ) is "homogeneous", the pres- 
ence of a non-trivial interaction singles out the loci (xi, . . . , xjv) G where the 
interaction cannot be neglected. In the physically relevant case of interactions rapidly 
decaying at infinity, the "diagonal" subsets like 

:= {x = (xi, . . . ,x N ) : xj = x k }, 1 < j ^ k < N, 

and particularly the "principal diagonal" 

D := {x= (x,x,...,x) : x e % d } 7 

along with their neighborhoods, are the zones where the inter-particle interaction is 
concentrated. Respectively, far away from all "diagonals" the interaction is cither 
zero (if the range of interaction is finite) or weak, thus negligible in the framework 
of some perturbative analysis. It turns out that a tedious geometric combinatorial 
analysis of the "diagonal" zones is not necessary for establishing qualitative results 
of the multi-particle Anderson localization theory. As was illustrated in [10], a much 
simpler classification of balls B^(u) suffices here: 

• balls close to the " principal diagonal" B ; 

• balls sufficiently distant from D. 

In the space Z N of indistinguishable configurations x of N particles, occupying 
pairwisc distinct positions xi , . . . , Xn € Z, we can no longer refer directly to the 
"diagonal" B. However, it is not the proximity to this set which is important per se 
for our analysis, but an easily verifiable fact that, for a point x = (xi, . . . ,xjv) € Z 

dist(x,B) < diam{xi, . . . ,x N } < 2dist(x,B). 

The quantity diamx = diam{xi, . . . ,xjv} is non-ambiguously defined for a configura- 
tion x = {xi, . . . , xn} of indistinguishable particles. 

Definition 4.1. An A^-particle ball B^(u) is said partially interactive (PI, in short) 
if diamu > A^L, An = 4iV, and fully interactive (FI), otherwise. 

The local Hamiltonian Hi , s with a finite-range interaction in a PI ball can is 

U 1,(11) ° 

algebraically decomposable in the following way: 

with ri + n" = N, u = (u', u"), u' E Z dn ' , u" G Z dn " , while for a fully interactive ball 
one cannot guarantee such a property. (In the case of a rapidly decaying interaction, 
the RHS is close in norm to the LHS.) 

Lemma 4.2. (A) //Bi(x) is PI, then there exist non-empty complementary subsets 

J,J C C [1,N] such that p(n J B L (x),n j7 cB L (x)) > 2L. 
(B) 7/Bi(x), B L (y) are two FI balls with p(x,y) > C N L, C N = 11N(> 2A N +3), 

then p(nB L (x),nB L (y)) > 2L. 

Proof. (B) Consider the set nB 2L (x) = Uj =1 B 2L (xj). If diamnx > A N L, this 
union cannot be connected, for otherwise it would have diameter < N ■ 4L, but 
diamIIB2L(x) > diamllx > A^L = ANL, Therefore, IIB i i+a(x) can be decomposed 
into a disjoint union of two non-empty clusters, 



J jeJ Bai(x i ))ll(Ui£j'» B 2L (xi)), 
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so for all j € J, i € J c one has 

p(B L (x 2 ),B L (^-)) > 2L, 

yielding the assertion (B). 

(C) If jo(x, y) > CnL, then for some j , p{xj ,yj ) > CnL. Since both balls are FI, 
for all j £ [1, n] 

p(xj,Xj ) < diam Hx < A N L, p{yj,Vj a ) < diam Lly < A N L, 
so by triangle inequality, for all i, j G 

p(xi,Vj) > p(x jo ,yjJ - p{%i,%j ) - p{Vj,yjJ > C N L~2A N L > 3L, 

yielding 

p(UB L (x),UB L (y)) >3L-2L = L. 

□ 

Note that in the case where a decomposition of the form (|4.6I) is possible, it may be 
not unique. For example, with N = 3 and an interaction of range r*o = 1, a particle 
configuration x = {0, 10, 20} can be decomposed into two non-interacting subsystems in 
three different ways. For our purposes, it will suffice to assume that one decomposition 
is canonically associated with every partially interactive ball. 

Now we state a result on localization in PI balls for Hamiltonians with a short-range 
interaction. The reader may want to compare the elementary proof of Lemma 
below with that of Lemma 3 in |10| which is more involved. 



Lemma 4.3. Assume that the interaction U has finite range r^ 6 [0, +oo). Consider 
an PI N -particle ball with canonical decomposition B^^ (u) = B^™ ' (uj) x B^™ '(ujc). 
If the balls B^" ^ (u') and B^™ ^ (u") are m—Coc, then B^(u) is also m—£c 



,oc. 



Proof. Since the interaction between subsystems labeled by J and by its complement 
J c vanishes, the operator H admits the decomposition 

H BMU ) = H <W) is + !<«') ®H <W) , (4.7) 

so that the eigenfunctions *f?j of H-B Lk ( u ) can be chosen in the form Sl/j = ® ipj„ , 
where tp^ are eigenfunctions of H ( „'). and, respectively, xpj,, are eigenfunctions of 

H („») . Now the required decay bounds for the spectral projections 

| VjMVjiy) | = | Vf ( X j) Vf ,(^)^r(yj)^f(y^) I (4.8) 

follow directly from the respective bounds on ip and xp, which are guaranteed by the 
m—Coc hypotheses on the projection balls Ty£'(u.j) and B^" \\ijc). □ 
Given numbers p, b > 0, introduce the following sequence: 

P(n,k) = P(n,k,p) = 2 N - n p(l + b) k , n = 1, . . . ,N; k > 0. (4.9) 



Lemma 4.4. Suppose that for all n G [l,iV — 1] and for a given k > the following 

-tide ball B^(u): 

| B^(u) is m-NCoc} < L~ k 



bound holds for any n-particle ball B^(u) 



-P(n,fc) 
J k 
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Then for any N -particle PI ball B^'(u) one has 

p{b£J(u) is m-MCoc) < 2L~ k 2P{N - k) < ±L* P(JV ' fc) . 

Proof. It follows from Lemma [4.31 that a ball (u) admitting a decomposition of 
the form (|4.7[) is m~M Hoc only if at least one of the projection balls Bj^-^uj-), 
B^" \ujc) is m—N Hoc. Therefore, with Lq large enough, we obtain 

p{b^ } (u) is m-JV£oc} 

<P^B { l l p(uja) is m-HCoc} +p{B^"" ) (u 7 c) is m-MCoc} 

< L- F{n '- k) + L - k P{n "> k) < 2L- 2P{N - k) < \L- k P{N > k) . □ 

4.4. Fully interactive iV-particle balls. Owing to Lemma 14.41 it suffices now to 
establish localization bounds for TV-particle FI balls Bj, fe+1 (x), assuming, if necessary, 
similar bounds 

• for n-particle balls of any radius L}-i, k' > 0, with n < N, 

• for iV-particle balls of any radius L k ' with k' < k. 

Lemma 4.5 (Main inductive lemma). Suppose that for a given k > and all i £ [0, k] 
the following bound holds for any N -particle ball B^°(u): 

P { b£°(u) is m-NCoc } < L; P{N ' l ' p \ (4.10) 

Then for any N-particle FI ball B^™' +i (x) one has 

p{bW i( u) ism-NCoc} <L- k P ^ k+ ^\ (4.11) 

provided that 

p>^jNd, 0<3&<min|^«i-Md ! V2-lj. (4.12) 
Furthermore, for any pair of distant balls Bj^(u), Bj^(v), one has 

P{]££l: B[ N J(u),B { £\v) are{E,m)-S} <L- p{N ' k > p) . (4.13) 

Proof. Set 

N k+1 = {B^ t (u) is m~N£oc}, 
S {2) = {3E3 distant (E,m)-S balls bJ»(u), b£?(v) C B^fw)}, 
U k 2} = {3E3 distant E-PR balls B^J (u),B { "J (v) C B^(w)}. 
By virtue of Lemma T3.61 we have A4+1 C TZ^ U (sj?^ \ , and by Theorem 14. 1[ 

P | TZ k 2) | < e" L fc , so that it remains to assess the probability P { S {2) \ TZ k 2) |. Fix 
points u,v £ B^^(w) and introduce the event (figuring in assertion (|4.13|) ) 

S {2) (u,v) = {3E : B^ } (u) and b£°(v) arc (E,m)-S }. 

Within the event <sf } (u, v)\^ 2) , either B^°(u) or B^(v) must be i?-CNR ; w.l.o.g., 
assume that B^p(u) is _E-CNR. Since it is (E, m)-S, Lemma [33] implies that it must 
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contain a pair of m-HCoc balls B^Jy'), B^fr") with d (y'>y") > h L l-i- Coil ~ 
sidcr two possible situations. 

(1) Both B^jy') and B^jy") are El. 

The supports nBi fe l (y'), nB£ fc x (y") of these of fully interactive balls have diam- 
eter smaller than \L)^ S (cf. Lemma |4~21 assertion (A)), so for Lq large enough, 

dist(nB Lfc (u),ILB ifc (v)) > Ll+ 5 - 2 ■ \L\ +S > \L\ +S . 
Using the inductive assumption (|4.10p and the mixing condition W2, we can write, 

1 /a 2 

with Lk-i = Lp!+i and Lq large enough: 

P{ B^»(y') and B^»(y") are m-M Hoc } < 2^ + f (1+ " )fc " 

The number of pairs y',y" is bounded by |B^ V -'(w)| 2 /2 < ^-^—L^ Nd , so in this case, 

p{ sf \^)}<^ t -S<«>'-'«« (4 . 14) 

A straightforward calculation shows that for Lq large enough, the RHS is bounded by 
jL k ^ N ' k+1 ' p \ under the conditions (|4.12[) . Indeed, we need the inequality 

^2 N ~ N p(l + bf' 1 - 2Nd > 2 N ~ N p(l + b) k+1 (4.15) 

or, equivalently, 

(1 + b) 2 < |§ - • 
Observe that for b G (0,1), (1 + b) 2 < 1 + 36, and ^ > p[ ™ b ) k -i ■ Therefore, if 
p > i^jNd and 36 < - then |g35| is also satisfied, so that 

p{sf }<p{^ 2) }+p{^ 2) \^ 2) } 

^ e + 1 L k+l < 2 L k+l 

yielding the assertion (|4.11[) in the case (1). 

(2) Either B^(y') or B^J^y") is PI. Then by Lemma Ol 

P{ Bi^fy') and B^Jy") is m-NCoe \ 

< min (p{ B^^y') is m-NCoc X ,P{ B^^y") is m-MCoc X) 

T -2 P (i+b) k - 1 T - P (i+b) k+1 

provided that (1 + b) 2 < 2, i.e., b < \[2— 1. This completes the proof. □ 
Applying Lemma 14. II and Lemma 14.41 we come by induction to the following result. 

Theorem 4.2 (Localization at any scale). Assume that \g\ is large enough, so that 
(|4TT ]l -P~2 ]) hold true with p > ^sNd. Then for any k > and any u e 7L Nd 

p{bW(u) m-AA£ 0C } < \L- k P ^ = (4.17) 

with b > satisfying (|4.12j) . In other words, with probability > 1 — ■|L fc p ^ 1+t '' , /or 
every eigenfunction tyj of operator H^n)^ and /or a/^ x,y G B^ (u) suc/i £/ia£ 



(4.16) 
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1 + iS + P 

|| x — y|| > L k ° one has 

This marks the end of the direct scaling analysis of localization of cigenfunctions in 
arbitrarily large finite balls. In Section El we will derive from the results of the scaling 
analysis uniform bounds on EF correlators in arbitrarily large finite balls. 

5. Scaling analysis: Adaptation to infinite-range interactions 

In this section, we adapt a certain number of definitions and statements to rapidly 
decaying, infinite-range interactions satisfying the assumption Ul. 

(51) Given a number S £ (0, 1/14), figuring in Ul, set 

q = 26, a = 1 + 46, t = 5/2, 

and observe that a 2 < 2, g — 5 > r. 

(52) replace the decay exponent 7(771, L) = m(l + L~ T ), figuring in the definition 13.31 
(cf. Eqn J232J0, by 

7 (m, L, 77) := m(l + L- T ) N - n+1 . 

(53) Modify the definition of distant balls as follows: balls B^(x), B^(y) are called 
distant iff p(x,y) > C N L 1+S , C N = 47V. 

(54) Modify the definition of PI and FI balls: a ball B L (x) is called PI if diam(x) = 
diam(xi, . . . , xm) > L 1+s , and FI, otherwise. 

(55) Check that if a ball (x) is PI, then it admits a decomposition 

B W(x)=BfV)xBfV) 

with dist(x',x") > 2L 1+S . 

(56) For every positive integer R, introduce the following finite-range approximations 
U_r, R > r , of a given interaction U generated by 2-body interaction potentials 
U^(r): 

Ug\r) = l r < R UW(r). 

(57) In order to measure the interaction between subconfigurations x', x" of an A^- 
particle configuration x £ , introduce the quantity 

e(R):= sup |U(x) - U(x') - U(x")|, x = (x',x"), 

d(x',x")>-R 

and check that for Rk = CnL^ +S , one has 

e(R k )<e-^ <\e- L ". 

The most important technical modification is required for the Lemma 14.31 the proof 
of which was very elementary, due to the use of eigenfunctions. 

Next, introduce the following finite-range approximations \Jr, R > tq, of a given 
interaction U generated by 2-body interaction potentials (r) : 

4 2) w = wc/ (2) m- 

Then it is readily seen that for any configuration x = (x',x") with dist(x',x") > R, 
one has U,r(x) = U^(x') + Unfa"). Furthermore, for any lattice subset of the form 
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B = Bl'(x') x B L »(x") with dist(B I/ (x'), Byr(x")) > R, the Hamiltonian 
admits the algebraic decomposition 

= Hg >(x>) <g> + l {n "> ®Hg //(x , f) , x' g Z dn ',x" g Z d "". 

Next, we fix an arbitrary m > 1 and modify the definition of a PI ball: B^ (x) 
is called PI iff diam(x) > Cn = 4l ^ t jv , and, therefore, it admits a decomposition 
Bf°(x) = Bz/(x') x B L »(x") with dist(x',x") > X = so that, for E fc = 

T^Lfc, one has 

c(iZfc) <e- 2mL " < ±e- L K 

Lemma 5.1. Assume that the interaction U satisfies the condition Ul. Fix an energy 
E. Let an N -particle ball B^ (u) = B^^Uj-) x B^™ \ujc), and a sample of the 
random potential V(-;lo) be such that 



(a) dist(x',x") > R k := %ftL k ; 

(b) B^'V) andB^'V) arem-Coc; 

(c) b£°(u) is £-CNR. 

Then the ball B^°(u) is (£?,m)-NS. 

Proof. Set i? fc = f^Lfe and define U (flfc) and H^ fc ) as above. Denote 

'J,:!' V = H^ )(u) , G(E) = (H B (jv )(u) - E)-\ G^(E) = (hJJ^ - 

Operator H^'^ admits the decomposition 



HS ,Ar) = H b<;')(uo ® + l(n ' } ® H b<;')(x»)' u = ( u '' u ") G ^ x Z * 1 " 

thus its eigenvalues are the sums E a ^ = \ a + fib, where {A a } = E(Hg i ( u /)) is the 

spectrum of H Bifc(u /) and, respectively, {fib} = S(H Bi) _(u"))- Eigenvectors of ~H.^ 
can be chosen in the form \I> a 5 = <p a ® if) b where {ip a } are eigenvectors of H 

and {rp b } are eigenvectors of H ( „"j . „.. Note that for each eigenvalue E a t, = X a + fib, 

i.e., for each pair (A a ,/it,), the non-resonance assumption \E — (A a + fib)\ > e Lk reads 
as \(E — X a ) — fib)\ > e L fc and also as \(E — fib) — A a )| > e L *=. Therefore, we can write 



= £ P' a (u', y') Gg-^j (u", y"; £ - A a ) (5.1) 
= J2 n (u", y") G<£> (u', y'; £7 - W ), 



where the resolvents G^ '^ u ,^(E — fib) and G^ ^ u „^(E — A ) arc non-resonant: 



\^!w&-l*)\\<<&, IIGb^V^B-A^II < 
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For any y 6 <9~B£ fc (u), either p(u',y') = Lk, in which case we infer from (b) 

|G (/?fc) (u,y;£)| < |B Lfc (u")|e-''( m ' i ^ Ar - 1 ) Lfc+L ^ (5.2) 
or p(u",y") = Lk, and then we have, respectively, 

\G (R »\u,y; E)\ < \B Lk (u')| e -« m ' L - N - 1 '> L * +L Z . 
In either case, the LHS is bounded by 

cxv(-m{l + Ll T ) N -( N -^ +1 L k + Ll + Const In L fc ) < le"^™'^). 

Next, by the second resolvent identity, G = G (i?fc) - G (Hfc) (U - U^ fc ))G, and using 
the assumed E-NK property of the resolvents G, G^', wc conclude that 

|| G- G^\\ < ||U-U( Hfr '|| ||G( Rfr )|| ||G|| < e - 2mLk e 2L t < i e -7(m,i»,A0. 

Now the claim follows from the inequality 

VxjeB^H |G(x,y ;J B)-G(^)(x,y;£;)|<Ie-^^). Q 

Lemma 5.2. Let B^ fc+1 (x), Bi fc+1 (y) be two distant PI balls. Then, for Lq large 
enough, 

P { B ifc+1 (x) and B Lk+1 (y) are (E, m)-S } < \L^ N ' k+1) . (5.3) 

Proof. Let 

TZ {2) ={3E: B Lfc+1 (x) andB Lfc+1 (y) are S-R } 

= { B Lfc+1 (x) and B Lt+1 (y) arc (E,m)-S } 

Then P { } < P { TZ^ } + P { \ TZ^ } , and within the event S (2 ~> \ TZ^ one of 
the balls Bi t+1 (x), B£ fc+1 (y) must be £?-NR, no matter how E £ M. is chosen. 
Next, consider the canonical decomposition of the PI ball 

B ifc+1 (x)=B ifc+1 (x')xB ifc+1 (x") 

By Lemma ??, if both Bi t+1 (x') and Bi fc+1 (x") arc m—Coc and Bi t+1 (x) is (E,m)- 
S, then Bi fc+1 (x) is E-K, thus B/, fc+1 (y) is i?-NR. On the other hand, consider the 
canonical decomposition of the PI ball 

B ife+1 (y)=B ife+1 (y')xB ife+1 (y"). 

If both B Lfc+1 (y') andB Lfc+1 (y") are m-Loc and B Lfi+1 (y) is (E,m)-S, then B^^y) 
must be E-K. Let 

£< 4 > = { one of the balls B Lfc+1 (x'), B Lk+1 (x"), B Lk+1 (y'), B Lfc+1 (y") is m-Af£oc}. 

Since these four balls correspond to systems with < N— 1 particles, we can use induction 
in N and write 

P{£< 4 > } < 4L^ 1 (Ar - 1 ' fe+1) = 4L^f (JV ' fc+1) . 
Then, as we have noticed, \ C TZ^- 2 \ hence, for Lq large enough, 

P { } < P { } + P { £< 4) } < e-^S + 4L^f (JV < fe+1) < Ii-^ fc+1 ). 

□ 

The statement of Lemma 13.51 remains unchanged, but its proof requires a minor 
modification. 
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Proof. ( Lemma \3.5\ under the hypothesis Ul.) (A) By assumption, either Hi, (u) is 
S-CNR and contains no (E, m)-S ball of radius L k , or there is a point w € Bj, fc+1 (u) 
such that any ball B ifc (v) C B Lfc+1 (u) with p(w,v) > CVi^ 5 is (E,m)-NS. In the 
former case, such an exclusion is unnecessary, but in order to treat both situations with 
one argument, we can formally set w = u (or any other point). 
Fix points x, y with R := p(x, y) > L h +e . By triangle inequality, 

p(x, B {CN -i) Lk (w)) + p(y, B (Cjv _ 1)Lfc (w)) >R- (2C N - 2)L k +s . 
Assume first that 

r' := p(x, B (Cn _ 1)l i+s (w)) > L k + 1, 

r" :=p(y,B {CN _ 1)Ll+ s(w))>L k + l. 

All balls of radius L k both in B r /(x) and in B r »(y) are automatically (E,m)-~NS. 
Furthermore, r' + r" > R- 2(C N - l)Ll +s -2>R- 2C N Ll +s . 

Consider the set B = B r /(x) x B r »(y) and the function / : B — > C defined by 
/(x', x") = Gb L;j+1 (x', x"; E). Since E is not a pole of the resolvent Ge 2jfe+1 (u)( - ); ^ is 
well-defined (hence, bounded, on a finite set). By Lemma l3.3l / is (Lk, q)-subharmonic 
both in x' and in x", with q < e _7 ' m ^ Lfc '"'. Therefore, one can write, with the conven- 
tion — InO = +00, using Lemma 13.21 and setting for brevity J := M — n + 1: 



ln/(u,y) > - In 



R-2C N L 1 k + S - 



e k + 1 



— ix yi - i, k ) - ^ k+1 ,-5 4) 

mR 



(1 + I^D J (1 - L, 1 ) (1 - 3C^«) - % 



> m (1 + r ) J > 7 (m, L fe+ i)i?. 

If r' = (rcsp., r" = 0), the required bound follows from the subharmonicity of the 
function /(x',x") in x" (resp., in x'). 

(B) This assertion is proved in the same way as in Section l3~3l □ 

The rest of the scaling procedure presented in Section^ does not have to be modified 
and applies to infinite-range interactions. 

6. Strong dynamical localization in finite volumes 

6.1. Uniform bounds in finite volumes. We use here a finite- volume variant of 
a 'soft' argument proposed by Gcrminct and Klein in |21j ; working with finite cubes 
allows to avoid a functional-analytic complement concerning the weighted Hilbcrt- 
Schmidt norms of spectral projections of operators H(w) in the entire graph and 
replace it with a simple application of Bessel's inequality. 

Denote by 3B\(T) the set of all Borel functions (f> : M — > C with supp <f> C I and 

Halloo <!■ " 

Theorem 6.1. Fix an integer L £ N* and assume that the following bound holds for 
any pair of disjoint balls Bi(x), Bi(y): 

P{3E G I : B L (x) and B L (y) are (E,m)-S) < f(L). 
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Then for any xj£ Z™ with d(x, y) > 2L+1, any connected subset A D Bi(i)UBi(j/) 
and any Borel function <fi G 3$\ (7) 

E[(|l x |0(H A H)|l y )|] <CL d c- mL + f(L). (6.1) 

Proof. Fix points x, y G Z™ with d(x, y) > 2L + 1 and a finite connected graph 
A D Bi(x) U Bx,(y). The operator H\{u>) has a finite orthonormal eigenbasis 
with respective eigenvalues {A^}. Set S = <9Bl(x) u 9Bx(y) (recall: this is a set of 
pairs (u, u')). Suppose that for some u, for each i there is z 6 {x, y} such that B^(z;) 
is Aj,m)-NS; let {v,} = {x,y} \ {zj. Denote A* x , y (</>) = (| l x \<j>{H\(u))\ l y )|, with 
Hx.,y{4>) < 1- Then by the GRI for eigenfunctions, and by Bessel's inequality used at 
the last stage of derivation, 

Mx, y (0) < Moo E l*»(x)* t (y)| < E l*<(^)*i(v<)l 
< E l*i(v.)|e- roi E |*i(u)| 



A,e/ (u,u')e9B L ( Zj 

E E i 

AiG/ (u,u')SS 



< e_mL E E i*i(u)i(i*i(x)i + 1**^)1) 



< e"™ L |S| max E \ (l**(u)| 2 + |*,(x)| 2 + |^(y)| 2 ) 

A s e/ 

< e -mL\S\ ^ 2|| ^ ||2 + || ^ ||2 + || x ||2) = e -mL| S | . 3 

where |S| < CL d . Denote 7V L = { 3E G I : B L (x) and B L (y) arc (E 1 , rn)-S }, with 
P{ A/"l } < /(£), by assumption. Further, 

E[ Mx , y (<£)] =E[l^Mx )y W]+E[l^/Xx,y(^)] <f(L) + 2CL d e~ mL . □ 

6.2. Strong dynamical localization in the entire graph Z™ . Here we follow the 
same path as in earlier works by Aizenman et al. (cf., e.g., [2J, [5]). 

Theorem 6.2. Consider the Hamiltonian H(w) of the form p. ID with random poten- 
tial satisfying the assumptions W1-W3 and the interaction potential satisfying one 
of the assumptions UO, Ul. Fix an interval I C K. There is go < +oo suc/i i/iai i/ 
M > go, then for all x, y G Z d , x ^ y, and with the same c, a > as in (??). 



3 



sup 

0lloo<l 



(hh)p / (h( w ))| i y ; 



< Constc- alnl+Cd ( x < y >. (6.2) 



Proof. For any ball B and any points x, y G B introduce a spectral measure /Ltg^ 
uniquely defined by 

/ 0(A)^B y JA) = (lxi^ffBHJnj^BMJliy), 

where is an arbitrary bounded Borel (or continuous) function, and similar spectral 
measures /x x , ,y for the operator H(cj) on the entire graph Z™ . If {BL fc } is a growing 
sequence of balls, then ^g' y u converge vaguely to fi*' y as k — > oo. Note that this 
fact remains true in a much more general context of unbounded operators: by a well- 
known result (cf., e.g., [22), for the strong resolvent convergence of operators H„ — ► H, 
with a common core T>, it suffices that H„</> — > H0 for any clement <p G 2?; in turn, 
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this implies the vague convergence of the spectral measures ((p,(f)(ll n )ij>). Usually, an 
appropriately chosen subspace of compactly supported functions can serve as a core, 
and on such functions operators H„ converge by stabilization. 

So, by virtue of Fatou lemma on convergent measures, for any measurable set £ C R 
and any growing sequence of balls Bi fc , 

E[|^|(£)]<liminf e[|^ \(£) 

Therefore, the uniform bounds in finite volumes Bi fc , established in Theorem ??, imply 
the dynamical localization on the entire lattice. □ 

Taking functions 4>t '■ A i— > e ttx , iel, one infers from Theorem l6.2l a more traditional 
form of dynamical localization: 

Theorem 6.3. Under the assumptions of Theorem \6.2\ there exist a,c> such that 
for any finite subset K C If 1 and any finite interval I C R 



E 



sup 

teR 



■ e -im*>) Pl (H(w))i 



K 



< oo. (6.3) 



Taking into account RAGE theorem(sjl, Theorem 16.31 implies the spectral localiza- 
tion: with probability one, the spectrum of H(w) is pure point, so that spectrally a.e. 
cigcnfunction \& of H(w) is square-summable. However, (|6.3[) does not imply directly 
an exponential decay of 

6.3. Exponential decay of eigenfunctions on the entire graph Z™ . 

Theorem 6.4. For P-a.e. lj £ f2 every normalized eigenfunction \& of operator H(w) 
satisfies the following bound: for some R(lo), x(oj) and all y with ||y|| > R(ui) 

|*(y)|<e-H|y|l. (6.4) 

Proof. Pick an arbitrary vertex z € Z™ . By Borel-Cantelli lemma, there is a subset 
Q! C f2 with P { f2' } = 1 such that for any w G f2' and some ko(uj), all fc > fco and any 
£ £ 1 there is no pair of Ll+ 5 -distant (E, m)-S balls B Lfc (x), B L& (y) C B Lfc+2 (z). 
Fix cj e Q'. 

Let * be a normalized eigenfunction of H(cj) with eigenvalue A. Since 11*112 < 1 
implies H^Hoo < 1, there is a point x such that \\^\\2 = |*(x)|. If x n € Bj, i ._ 1 (0), 
then Bi fc (0) must be (A, m)— S, otherwise the (A, m)-NS property would lead to a 
contradiction: 

ll*lloo = |*(i)|<e-^||*|| 00 <||*|| 00 . 
Thus any ball B Lfc (y) C B Lfc+2 (z) with p(y,z) e [Lfc+i, Lfc+2) (hence, with p(y, z) > 
L£ > 2iJ, + ' 5 ) is (A, m)— NS. Note that the function x i-> |*„(x)| is (L kl q)-subharmonic 
in B/j(y), with i? = p(y, z) - 2L^+' 5 - 1 > \L 1 k +2S and g = e -i(™,L k )L K Rccall that 
we set t = (5/4; cf. (|3.19[) . Now Lemma |3~T1 implies, for Lk large enough, 

> - (1 + L D (l " f^) > - (1 + * V) > - 
yielding the assertion (|6.4[) . □ 



''See, e.g., the original papers 131 1161 and their discussion in 1121 
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7. Fermionic Hamiltonians on more general graphs 

The reduction to a standard lattice Laplacian on a subset {(xi , . . . , xn) : x\ < ■ ■ ■ < 
xn} with Dirichlct boundary conditions is no longer possible for particle systems on 
lattices 7L d with d > 1. Instead, one has to work with a symmetric power of the lattice, 
considered as a graph. So it seams reasonable to consider a fairly general, countable 
connected graph Z satisfying the condition of polynomial growth of balls: 

Vxe2 VL > 1 |Bi(a:)| < C d L d , 

where Hl(x) — {y £ Z : dz(x,y) < L}. In particular, this gives a uniform bound on 
the coordination numbers, nz(x) < Cd (of course, this bound may be non-optimal). 

An orthonormal basis in the Hilbert space of square-summable antisymmetric func- 
tions : Z N — > C is formed by the functions 

1 x -, N 

* a = -7= > (g) l a , , a = {ai, ...jOjv}, #{ai,...,ojv} = N. 
VN\ f£ 3=1 * <J) 

Further, define the graph {ZU,SU) as follows: the vertex set is 

Z> = {a = {01, . . . , a N } : Oj G Z, #{01, . . . , 0^} = N}. 

Two vertices a, b form an edge iff 

• the symmetric difference a0b has cardinality 1, i.e. a = {ai, c%, . . . , cjv}, 
b = {61, c 2 , . . • ,cjv}, with #{ai,6i,c 2 , . . . ,cat} = N + 1, and 

• d 2 (ai,6i) = 1. 

Next, define on the max-distance, 

p(x,y)= min max . d 2 (x ff -i (j) , 

and introduce the balls Bj,(x) relative to the distance p(- , •). 

Now one can define the fermionic (negative) Laplacian (—A) on Z™ and random 
Hamiltonians H(w) = — A + gV(x; w) + U(x), where 

U(x)=5^C7'(d(a!i > Xj)), C/:N^M, 

and the external random potential energy 

V(x;o;) = V{x 1 ;oj) + --- + V{x N ]uj), 

is generated by a random field : Z x SI — >• M. 

The method presented in Sections [3]-[5] applies to strongly disordered random Hamil- 
tonians H(w) = —A + jV(x;w) + U(x) describing fermionic systems on connected 
graphs Z with polynomial growth of balls; indeed, one can see that we did not use 
particular properties of the one-dimensional lattice Z = Z. 

Strongly disordered bosonic systems can be treated in a similar way; the only modifi- 
cation required here concerns the explicit form of the matrix elements of the Laplacian, 
which remains a second-order finite-difference operator. 
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Appendix. Proofs of auxiliary statements 

7.1. Proof of Lemma 14.11 We start with the second assertion. 
The random potential energy V(<j) reads as follows: 

V(x 1 ,...,x N ;u) =^j\yV{y\u). (7.1) 

2/Sx 

Therefore, if p(x,y) ^ 0, then there exists a point w € IIBi(u) such that n m (x) ^ 
n w (y). As a result, 

V(x;w) - V(y;w) - (n m x) - n^y^V^u;) + £ CyV(v,u), (7.2) 

where the explicit form of the integer coefhcients c v = n„(x) — n t ,(y) is irrelevant for our 
argument: it suffices to know that the sum in the RHS of Eqn. (|7.2p is measurable with 
respect to the sigma-algebra generated by the random variables {V(v; •), v ^ w}, 
while n m (x) — n w (y) =: c w ^ 0, \c w \ > 1. Therefore, 

P{|.gV(x; W )-.gV(y; W )| < S }=E[P{|V(x; W )-V(y; W )| < Ig^s^}] 

= E[V{\c u ,V(w;u)+((lj)\ < Igp 1 !^}] 
with some S^^-measurable random variable 

= E [P { c w V(w; u) e [CM - Iffl" 1 ^ C(w) + I } ] • 

Since |cu,[ > 1 and the interval [£(u>) — s, + s] has length 2s, we conclude that 
P{|V(x;w)-V(y;w)| < {g^s} < sup (F v , w (a + 2s) - F v , w (a)). 

Since Fy tW is continuous, by assumption Wl, the latter quantity vanishes as 4- 0. 

Therefore, 

P{3x,yeB L (u) : x^y, | fl V(x) - <?V(y)| < a } —> 0, 

so, with arbitrarily high probability, the spectrum of the diagonal operator V(w) in 
Bi (u) admits a positive uniform lower bound s > on spectral spacings (differences 
between the eigenvalues). By taking \g\ large enough, all spacings for operator g~V(uj) 
can be made arbitrarily large. Eigenvectors of a continuous hnitc-dimcnsional operator 
family A{t) with simple spectrum at t = to are continuous in a neighborhood of to- To 
prove the second assertion, it suffices to apply this fact to the family A{t) = "W—g~HA., 

te [0,1]. 

The proof of the first assertion is even simpler. Using again the representation (|7.1I) . 
we see that for each x the value of the potential energy V(x;w) + U(x) is a linear 
combination (with integer coefficients) of random variables with continuous probability 
distribution obeying Wl. Arguing as above, one can see that, for any E <G R, s > 
and |<7| large enough, with probability arbitrarily close to 1, dist(£ l , S(H Bl/0 ( u ))) > s - 
By Combes-Thomas estimate, this implies exponential decay of the Green functions, 

|G Bio( u)(x,y;£)| < e - m ^ d ^ 

with m(s) — > oo as s — > oo. Since the graph distance d(- , ■) on dominates the 
max-distance p(- , ■), the claim follows. □ 
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7.2. Proof of Theorem 14. 11 We prove a result slightly stronger than required for the 
spectral analysis of fermionic operators H(w): an EVC bound is established for opera- 
tors in the Hilbert space of quantum states of distinguishable particles, not just in the 
subspace of anti-symmetric functions. This simplifies geometrical arguments and nota- 
tions. So, in this subsection, we denote, for each given n > 1, x = (xi, . . . , x n ), B^(x) = 
{y G (Z d )" : |y — x| < L}; A stands for the nearest-neighbor lattice Laplacian in (Z d ) n , 
and H B[ ( X ) = — A Bi ( x ) + c/V + U. We denote by ps stands for the symmetrized dis- 
tance: /Os( x i y) = m i n 7res„ |tt( x ) — y| where the elements n of the symmetric group 
act by permutations of the coordinates: tt(x±, . . . ,x n ) — (a^-im, . . . 

Given an Af-particlc configuration x and a proper index subset J C denote 
by Ilyx the set {x j} j £ J} C Z d . 

Definition 7.1. A cube B^(x) is weakly separable from B^(y) if there exists a par- 
allelepiped Q C Z d in the 1-particlc configuration space, of diameter R < 2NL, and 
subsets J\,J 2 C [1,N] such that \ > \J 2 \ (possibly, with J 2 = ) an d 

n Jl B L (x)un j2 B i ( y ) eg, 
n^.B L (y)ng = 0. 

A pair of balls (B^(x), B^(y)) is weakly separable if at least one of the balls is weakly 
separable from the other. 

Lemma 7.1 (Cf. Lemma 2.3 from [7]). Cubes Bi(x),Bi(y) with d(x, y) > ANL are 
weakly separable. 

See the proof in [7J . 

Lemma 7.2. Let V : 1 d x 51 — ► R be a random field satisfying the condition (W3). 
Let x, y G Jj Nd be two configurations such that the balls B^(x), B^(y) are weakly 
separable. Then for any s > the following bound holds for the spectra of operators 
Hb £ ,(x)> Hb „( y ) with arbitrary L',L" < L: 

P{dist(a(H Bl ,(x)) ! (T(HB i „( y ))) <s}< h L (2s). 
with h^ defined in (14.41) . 

Proof. Let Q be a parallelepiped which satisfies the conditions (|7.3|) . 

n Jl B L (x)un j2 B L ( y ) eg, 

n 7 =B L (y)nQ = 0, 

for some J\,3 2 C [l,iV] with \ Ji\ = n\ > n 2 = \Ji\- In terms of the sample mean £ = 
£q and the fluctuations {rj x , x € Q} defined in subsection 12.51 operators Hb ,( x )(w), 
H Bi „( y )(w) read as follows: 

H Bj ., (x) (w) = ni£(w) 1 +A(w), H Bj .„(y)(w) - n 2 f (w) 1 +B(w) (7.4) 

where operators A(w) and B(w) are SxQ-mcasurable. Let 

{EZ,..., E> M ,}, M'= |B L ,(x)|, 
{^',...,^„},M" = |B^(y)| 

be the sets of eigenvalues of H B ( x ) and of H B[ „( y )) counted with multiplicities. 
Owing to Eqn (|7.4p . these eigenvalues can be represented as follows: 

E' 3 {lo) = mtiw) + Efa), E'!{u) = n 2 £(Lo) + E%(u), 
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where the random variables E'^iui) and E"(uj) are SV^-measurable. Therefore, 

E[{u) - £» = <H + (£» - £»), 
with n := n\ — n-i > 1, owing to our assumption. Further, we can write 

P { dist (cr(H Bi , (x) ) , <7(H Bi „ (y) )) < a } = P { 3 i, j : \E[ - E'j | < a } 



A/' A/" 



A/' M" 



<^^P{|^-4'|< S } = ^^E[P{|^-4'|< S |^, Q }] 
i=l j=l i=l j=l 

For all i and j, we have 

P{|^-^ / |<s|ffv I g}=P{|(n 1 -n 2 )f + ^-^ / |<8|?v I o} 



E'!-E[ 



E" - E' 



< ^ L (2n- 1 s|5y,Q) < v l {1s\$ v ,q)- 

Therefore, 

P{dist( CT (H Bt „(x)MHB t „(y))) <«} <M'M"v L (\n 1 -n 2 \- 1 8) 

< \B L „(x)\.\B L ,,(y)\v L (2a). 
Now Theorem 14.11 follows from Lemma 17.21 combined with Lemma 17.11 



(7.5) 



□ 
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